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ABSTRACT 

We present the results from the Suzaku X-ray observations of five 
fiat-spectrum radio quasars (FSRQs), namely PKS 0208-512, Q 0827+243, 
PKS 1127-145, PKS 1510-089 and 3C 454.3. All these sources were addition- 
ally monitored simultaneously or quasi-simultaneously by the Fermi satellite in 
gamma-rays and the Swift UVOT in the UV and optical bands, respectively We 
constructed their broad-band spectra covering the frequency range from 10^^ Hz 
up to 10^^ Hz, and those reveal the nature of high-energy emission of luminous 
blazars in their low-activity states. The analyzed X-ray spectra are well fitted by 
a power-law model with photoelectric absorption. In the case of PKS 0208— 512, 
PKS 1127—145, and 3C 454.3, the X-ray continuum showed indication of hard- 



^^Institute of Space and Astronautical Science, JAXA, 3-1-1 Yoshinodai, Sagamihara, Kanagawa 229-8510, 
Japan 



^•^Corresponding author, rsato@astro.isas.jaxa.jp 

^^Institut fiir Astro- und Teilchenphysik and Institut fiir Theoretische Physik, Leopold-Franzens- 
Universitat Innsbruck, A-6020 Innsbruck, Austria 

^^Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA 

^^Department of Physics, University of Washington, Seattle, WA 98195-1560, USA 

^"^Santa Cruz Institute for Particle Physics, Department of Physics and Department of Astronomy and 
Astrophysics, University of California at Santa Cruz, Santa Cruz, CA 95064, USA 

s^NYCB Real-Time Computing Inc., Lattingtown, NY 11560-1025, USA 

^® Astronomical Observatory, Jagiellonian University, 30-244 Krakow, Poland 

^''Department of Chemistry and Physics, Purdue University Calumet, Hammond, IN 46323-2094, USA 

^^Institut fiir Theoretische Physik and Astrophysik, Universitat Wiirzburg, D-97074 Wiirzburg, Germany 

^^Consorzio Interuniversitario per la Fisica Spaziale (CIFS), 1-10133 Torino, Italy 

60JNTEGRAL Science Data Centre, CH-1290 Versoix, Switzerland 

^^INAF, Osservatorio Astronomico di Torino, 1-10025 Pino Torinese (TO), Italy 

^^Dipartimento di Fisica, Universita di Roma "Tor Vergata" , 1-00133 Roma, Italy 

^■^School of Pure and Applied Natural Sciences, University of Kalmar, SE-391 82 Kalmar, Sweden 

®'*Nicolaus Copernicus Astronomical Center, Bartycka 18, 00-716, Warsaw, Poland 

^^MuUard Space Science Laboratory/UCL. Holmbury St Mary, Dorking, Surrey RH5 6NT 

®^ Osservatorio Astronomico di Brena, via Brera, 28, Milano 1-20121, Italy 



-5 - 



ening at low-energies. Moreover, when compared with the previous X-ray ob- 
servations, we see a significantly increasing contribution of low-energy photons 
to the total X-ray fluxes when the sources are getting fainter. The same be- 
havior can be noted in the Suzaku data alone. A likely explanation involves a 
variable, flat-spectrum component produced via inverse-Compton (IC) emission, 
plus an additional, possibly steady soft X-ray component prominent when the 
source gets fainter. This soft X-ray excess is represented either by a steep power- 
law (photon indices F ~ 3 — 5) or a blackbody-type emission with temperatures 
kT ~ 0.1 — 0.2keV. We model the broad-band spectra spectra of the five observed 
FSRQs using synchrotron self-Compton (SSC) and/or external-Compton radia- 
tion (ECR) models. Our modeling suggests that the difference between the low- 
and high-activity states in luminous blazars is due to the different total kinetic 
power of the jet, most likely related to varying bulk Lorentz factor of the outflow 
within the blazar emission zone. 

37 Subject headings: galaxies: active — quasars: jets — radiation mechanisms: 

38 non-thermal - X-rays: galaxies 

39 1. Introduction 

40 Observations with the EGRET instrument (30 MeV to 30 GeV; Thompson et al. 1993) 

41 on board the Compton Gamma-Ray Observatory (CGRO) have resulted in detection of 7- 

42 ray emission from a few hundred astrophysical sources, 66 of which were securely associated 

43 with active galactic nuclei (AGNs; e.g., Hartman et al. 1999). Most of the AGNs detected 

44 by EGRET show characteristics of the blazar class. Obscrvationally, this class include flat- 

45 spectrum radio quasars (FSRQs) and BL Lac objects. FSRQs have strong and broad optical 

46 emission lines, while the lines are weak or absent in BL Lacs. During the flrst three months of 

47 the Fermi Large Area Telescope's (LAT) all-sky-survey, 132 bright sources at high Galactic 

48 latitudes (|6| > 10°) were detected at a confldence level greater than 10 a (Abdo et al. 

49 2009a). As expected from the EGRET observations, a large fraction (106) of these sources 

50 have been associated with known AGNs (Abdo et al. 2009b). This includes two radio 

51 galaxies (CentaurusA and NGC 1275; Abdo et al. 2009c) and 104 blazars consisting of 58 

52 FSRQs, 42 BL Lac objects, and 4 blazars with unknown classification based on their Spectral 

53 Energy Distribution (SED). 

54 The radio-to-optical emission of luminous blazars of the FSRQ type is known to be pro- 

55 duced by the synchrotron radiation of relativistic electrons accelerated within the outflow, 

56 while the inverse Compton (IC) scattering of low-energy photons by the same relativistic 
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57 electrons is most likely responsible for the formation of the high energy X-ray-to-7-ray com- 

58 ponent. In addition, it is widely believed that the IC emission from FSRQs is dominated 

59 by the scattering of soft photons external to the jet (external Compton radiation, ECR). 

60 These photons, in turn, are produced by the accretion disk, and interact with the jet either 

61 directly or indirectly, after being scattered or reprocessed in the broad-line region (BLR) or 

62 a dusty torus (DT; see, e.g., Dermer & Schlickeiser 1993; Sikora et al. 1994; Blazejowski et 

63 al. 2000). Other sources of seed photons can also contribute to the observed IC radiation, 

64 and these are in particular jet synchrotron photons through the synchrotron self- Compton 

65 process (hereafter SSC; Maraschi et al. 1992; Sokolov & Marscher 2005). 

66 In this context, detailed X-ray studies offer a unique possibility for discriminating be- 

67 tween different proposed jet emission models, since those scenarios predict distinct compo- 

68 nents to be prominent in blazar spectra around keV photon energies. For example, in the 

69 soft X-ray range a break is expected in the ECR/BLR model, tracking the low-energy end 

70 of the electron energy distribution (Tavecchio et al. 2000; Sikora et al. 2009). Indeed, both 

71 the XMM-Newton and the Suzaku X-ray data of RBS 315 show "convex" spectra (Tavec- 

72 chio et al. 2007). Such a curvature, on the other hand, can be alternatively accounted for 

73 by an excess absorption below 1 keV over the Galactic value, or by an intrinsic curvature 

74 in the electron energy distribution. Furthermore, the situation can be more complex, with 

75 the simultaneous presence of yet additional components, such as the high-energy tail of the 

76 synchrotron continuum, SSC emission, or the narrow-band spectral feature originating from 

77 the "bulk Comptonization" of external UV (disk) radiation by cold electrons within the 

78 innermost parts of relativistic outflow (Begelman & Sikora 1987; Sikora & Madejski 2000; 

79 Moderski et al. 2004; Celotti et al. 2007). 

80 GhiseUini et al. (1998) have studied the spectral energy distribution of 51 EGRET- 

81 detected 7-ray loud blazars and have apphed the SSC-I-ECR model to the spectra of these 

82 sources. Although most of the broadband data collected by GhiseUini et al (1998) cor- 

83 responded to non-simultaneous measurements, those authors discovered clear trends and 

84 correlations among the physical quantities obtained from the model calculations. In par- 
as ticular, they found an evidence for a well-defined sequence such that the observed spectral 

86 properties of different blazar classes (BL Lacs and FSRQs) can be explained by an increasing 

87 contribution of an external radiation field towards cooling jet electrons (thus producing the 

88 high-energy emission) with the increasing jet power. As a result, while the SSC process alone 

89 may account for the entire high-energy emission of low-power sources (BL Lacs), a significant 

90 contribution from the ECR is needed to explain the observed spectra of high-power blazars 

91 (FSRQs). Meanwhile, when focusing on one particular object, Mukherjee et al. (1999) 

92 reported that they found a similar trend in the different spectral states of PKS 0528+134. 

93 They studied the sequence of flaring and low-flux states of the source and found that the 
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94 SSC mechanism plays a more important role when the source is in a low state, and the ECR 

95 mechanism is the dominant electron cooling mechanism when the source is in a high 7-ray 

96 state (see in this context also Sambruna et al. 1997). 

97 In order to understand the blazar phenomenon and the differences between BL Lacs and 

98 FSRQs, as well as the origin of spectral transitions in a particular object, one has to obtain 

99 truly simultaneous coverage across the entire spectrum, diuing both flaring and low-activity 

100 states. However, past 7-ray observations in low-activity states have been limited to only a few 

101 extremely luminous objects, such as PKS 0528-134 or 3C279. Only now, with the successful 

102 launch of the Fermi satellite and the excellent performance of the Suzaku instruments, do 

103 we have an opportunity to study high-energy spectra of blazars with substantially improved 

104 sensitivity, and therefore can probe the different states of the sources' activity. 

105 In this paper, we report the high-sensitivity, broadband Suzaku observations of five 

106 FSRQs, namely PKS 0208-512, Q 0827+243, PKS 1127-145, PKS 1510-089, and 3C 454.3, 

107 which were bright gamma-ray sources detected by EGRET. Additionally, all of these sources 

108 were monitored simultaneously or quasi-simultaneously by the Fermi LAT and Swift Ultravi- 

109 olet/Optical Telescope (UVOT; Roming et al. 2005). These broadband and high-sensitivity 
no observations allow us to reveal the characteristics of the high-energy IC continuum in the 

111 low-activity states of luminous blazars. The paper is organized as follows: in §2, we describe 

112 observation and data reduction in the X-ray {Suzaku)^ UV-optical (Swift UVOT) and 7-ray 

113 [Fermi LAT) domains. In §3, we present the broad-band analysis results. Finally, in §4 

114 we discuss the constraints on the jet parameters and speculate on the the origin of differ- 

115 ent activity states in luminous blazars. Throughout the paper we adopt the cosmological 

116 parameters Hq — 71kms'"^ Mpc~^, Om = 0.27, and Qa — 0.73. 



117 2. Observation and Data Reduction 

118 2.1. Suzaku 

119 Five FSRQs were observed by Suzaku (Mitsuda et al. 2007) for 40 ks each as one of 

120 the long-category projects between 2008 October and 2009 January. Table 1 summarizes the 

121 start time, end time, and the exposures for each observation. Suzaku carries four sets of X-ray 

122 telescopes (Scrlcmitsos et al. 2007), each with a focal-plane X-ray CCD camera (XIS, X-ray 

123 Imaging Spectrometer; Koyama et al. 2007) that is sensitive over the 0.3 — 12keV band, 

124 together with a non-imaging Hard X-ray Detector (HXD; Takahashi et al. 2007; Kokubun et 

125 al. 2007), which covers the 10 — 600 kcV energy band by utilizing Si PIN photo-diodes and 

126 GSO scintillation detectors. All of the sources were focused on the nominal center position 
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127 of the XIS detectors. 

128 For the XIS, we used data sets processed using the software of the Suzaku data process- 

129 ing pipehne (ver. 2.2.11.22). Reduction and analysis of the data were performed following 

130 the standard procedure using the HEADAS v6.5 software package. The screening was based 

131 on the following criteria: (1) only ASCA-grade 0,2,3,4,6 events were accumulated, while hot 

132 and flickering pixels were removed using the CLE ANSIS script, (2) the time interval after the 

133 passage of South Atlantic Anomaly was greater than 500 s, and (3) the object was at least 

134 5° and 20° above the rim of the Earth (ELV) during night and day, respectively In addition, 

135 we also selected the data with a cutoff rigidity (COR) larger than 6 GV. The XIS events were 

136 extracted from a circular region with a radius of 4.2' centered on the source peak, whereas 

137 the background was accumulated in an annulus with inner and outer radii of 5.4' and 7.3', 

138 respectively. We checked that the use of different source and background regions did not 

139 affect the analysis results. The response and auxiliary files were produced using the analysis 

140 tools XISRMFGEN and XISSIMARFGEN developed by the Suzaku team, which are included in 

141 the software package HEAsoft version 6.5. 

142 The HXD/PIN data (version 2.0) were processed with basically the same screening 

143 criteria as those for the XIS, except that we required ELV > 5° through night and day and 

144 COR> 8GV. The HXD/PIN instrumental background spectra were provided by the HXD 

145 team for each observation (Kokubun ct al. 2007; Fukazawa et al. 2006). Both the source 

146 and background spectra were made with identical good time intervals and the exposure 

147 was corrected for detector dcadtimc of 6.0 — 8.0%. We used the response files, version 

148 ae_hxd_pinhxdnom5_20080716.RSP, provided by the HXD team. In our analysis, the 

149 hard X-ray emission of PKS 1510—089 and 3C454.3 were detected in the energy range from 

150 12 keV to 40 keV and 50 keV, respectively. For other objects, the sources were not detected 

151 in the HXD/PIN data. We also note here that all of the objects, the sources were not 

152 detected in the HXD/GSO data. 

Table 1: Suzaku observation log of five FSRQs. 



Object 


z 


Start time 


Stop time 


XIS/HXD exposures 






(UT) 


(UT) 


(ks) 


0208-512 


1.003 


2008 Dec 14 07:33 


2008 Dec 15 11:30 


50.3/39.3 


0827+243 


0.939 


2008 Oct 27 05:11 


2008 Oct 28 08:04 


35.3/36.3 


1127-145 


1.187 


2008 Nov 29 18:10 


2008 Nov 30 22:51 


42.2/29.0 


1510-089 


0.361 


2009 Jan 27 04:32 


2009 Jan 28 05:25 


38.5/36.2 


3C 454.3 


0.859 


2008 Nov 22 09:19 


2008 Nov 23 16:31 


39.9/40.4 
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153 2.2. Swift 

Four analyzed FSRQs (PKS 0208-512, Q 0827+243, PKS 1510-089, and 3C 454.3) were 

155 observed with Swift between 2008 October and 2009 January, as part of Swift "target of 

156 opportunity" observations. We analyzed the data taken within or near the time of the 

157 Suzaku observations. For the case of PKS 1127— 145, however, the observations were made 

158 only once in 2007 March. We focused on analysis of the UVOT data, since Suzaku provides 

159 much better photon statistics in X-rays than Swift X-ray Telescope (XRT; Burrows et al. 

160 2005) and Burst Alert Telescope (BAT; Barthelmy et al. 2005), thanks to the long Suzaku 

161 exposures. We used the XRT data primarily for a consistency check regarding the spectral 

162 properties. Table [2] summarizes the start time, exposure time, and filters used for each 

163 observation. 

164 The UVOT observing mode commonly takes an exposure in each of the six optical and 

165 ultraviolet filters (f , b, u, uvwl, uvm2, and uvw2) per Swift pointing. The list of UVOT 

166 observations is given in Table [21 For the screening, reduction and analysis of the Swift data, 

167 we used standard procedures within the HEASoft v.6.5 software package with the calibration 

168 database updated as of 2009 February 28. For this analysis. Level 2 sky-corrected image data 

169 were used. Since all sources were relatively bright, the source aperture sizes were chosen to 

170 correspond to those used to determine the UVOT zero points: 5" for the optical and UV 

171 filters (Poole et al. 2008). The background was extracted from a nearby source-free circular 

172 region with 15" radius. All image data were corrected for coincidence loss. The observed 

173 magnitudes were converted into flux densities by the standard procedures (Poole et al. 2008). 

174 The XRT data were all taken in Photon Counting mode (PC mode; Hill et al. 2004). The 

175 data were reduced by the XRT data analysis task xrtpipline version 0.12.0. Photons were 

176 selected from the event file by xselect version 2.4. The auxiliary response file was created by 

177 the XRT task xrtmkarf and the standard response file swxpc0to12s6_20010101v01 l.RMF 

178 All spectra were analyzed in the 0.3—10.0 keV band using XSPEC version 11.3.2. 

179 2.3. Ferrm LAT 

180 During the first year of Fermi Large Area Telescope (LAT; Atwood et al. 2009) oper- 

181 ation, most of the telescope's time has been dedicated to "survey mode" observing, where 

182 Fermi points away from the Earth, and nominally rocks the spacecraft axis north and south 

183 from the orbital plane to enable monitoring of the entire sky every ~ 3 hours (or 2 orbits). 

184 We analyzed the LAT's observations of the five blazar regions using data collected during 

185 the first 4-5 months centered around Suzaku observations. Little variability indicated by the 
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Table 2: Swift observation log of five FSRQs. 



Object 


obsID 


Start time 


Exposure" 


Exposure^ 


Filter^ 






(UT) 


(ks) 


(ks) 




0208-512 


00035002024 


2008 Dec 14 15:25 


0.99 


0.94 


all 


0827+243 


00036375004 


2008 Dec 08 13:42 


1.72 


1.71 


u 


1127-145 


00036380001 


2007 Mar 24 00:32 


14.6 


14.2 


all 


1510-089 


00031173010 


2009 Jan 25 18:40 


3.46 


3.40 


u, wl, m2 


3C 454.3 


00035030030 


2008 Oct 26 20:28 


0.43 


0.40 


all 



"Swift XRT 
^Swift UVOT 

186 LAT lightcurves for the studied objects during this time implies that the constructed broad- 

187 band spectra, even though not exactly simultaneous, are representative for the low-activity 

188 states of all five blazars. 

189 The data used here comprise all scientific data obtained between 4 August and 19 

190 December 2008 for PKS 0208-512, Q 0827+243, PKS 1127-145 and 3C 454.3 (interval runs 

191 from Mission Elapsed Time (MET) 239557417 to 251345942), and 4 August 2008 and 30 

192 January 2009 for 1510-089 (MET 239557417 to 254966035), respectively . We have apphed 

193 the zenith angle cut to eliminate photons from the Earth's limb, at 105°. This is important 

194 in pointed mode observations, but also important for survey mode due to overshoots and sun 

195 avoidance maneuvers. In addition, we excluded the time intervals when the rocking angle 

196 was more than 43°. We use the "Diffuse" class events (Atwood et al. 2009), which, of all 

197 reconstructed events have the highest probability of being photons. 

198 In the analysis presented here, we set the lower energy bound to a value of 200 MeV, 

199 since the bin counts for photons with energies of ~ 100 MeV and lower are systematically 

200 lower than expected based on extrapolations of a reasonable functions. Science Tools version 

201 v9rl4 and IRFs (Instrumental Response Functions) P6_V3 were used. 



202 



3. Results 



203 3.1. Suzaku 

204 3.1.1. Temporal analysis 

205 Figure [T] shows the count rate variations of the five observed FSRQs. The summed XISs 

206 (XIS0,1,3) hght curves are shown separately in different energy bands: 0.5— 2keV [upper 

207 panel), 2— lOkeV {middle panel), and 0.5— lOkeV {bottom panel), respectively. Since the 
20B count rate variations of the HXD/PIN detector were less clear due to limited photon statistics 

209 and uncertainty of the modeling of the non-X-ray background, we only concentrate on the 

210 temporal variability of the XIS data below 10 keV. We evaluate the fractional variability by 

211 calculating the variability amplitude relative to the mean count rate corrected for effects 

212 of random errors (e.g., Edelson et al. 2002): Fvar = a/S^ — ct^^/x, where S"^ is the total 

213 variance of the light curve, cxerr^ is the mean error squared and x is the mean count rate. 

214 The variability amplitude in the XIS bands are F^ar = 0.036 ± 0.021 for 0208—512, F^ar = 

215 0.027 ± 0.010 for 1127-145, and F^ar = 0.025 ± 0.008 for 1510-089, respectively. 0827+243 

216 and 3C454.3 show only weak variability, which is not significant. 



217 3.1.2. Time- averaged spectral analysis 

21B In the following we report the analysis procedure and results for each object. The 

219 background-subtracted spectra were fitted using XSPEC ver.11.3.2. All errors are quoted at 

220 the 90% confidence level for the parameter of interest unless otherwise stated. All the fits in 

221 this paper are restricted to the energy ranges of 0.5 — lOkeV (X1S0,3: the FI chips), 0.3— 8keV 

222 (XlSl: the BI chip), 12-40 keV for PKS 1510-089, and 12-50 keV for 3C 454.3 (HXD/PIN). 

223 We fixed the relative normalization of the XISs and HXD/PIN at 1.13, which is carefully 

224 determined from the XIS calibration using nominal pointings of the Crab Nebula. Serlemitsos 

225 et al. (2007) reported that spectral normalizations are slightly different (a few percent) 

226 among the CCD sensors based on a contemporaneous fit of the Crab spectra. Therefore, 

227 we adjusted the normalization factor among the three XISs relative to XISO. The results of 
22B the spectral fits with a simple absorbed power-law model are summarized in TableE] (with 
229 Galactic absorption) and TablelH 
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Table 3: Results of the spectral fits to the Suzaku data using a power-law with Galactic 
absorption 



Object 




r 




0208-512 
0827+243 
1127-145 
1510-089 
3C 454.3 


3.08 (fixed) 
3.62 (fixed) 
3.83 (fixed) 
7.88 (fixed) 
7.24 (fixed) 


1.68±0.03 
1.46±0.04 
1.41±0.02 
1.37±0.01 
1.58±0.01 


1.37±0.06 
1.37±0.07 
3.45±0.08 
6.31±0.12 
16.7±0.2 



constant 

(XIS0,1,3,HXD/PIN) 



(dof) 



l,1.04±0.05,1.04±0.05,None 0.91 (250) 

l,0.90±0.05,1.04±0.06,None 0.84 (194) 

l,l-03±0.03,l-05±0.03,None 1.01 (331) 

1,L00±0.02,1.02±0.02,1.13 1.06 (407) 

1,L04±0.01,1.02±0.01,1.13 1.00 (1090) 



Errors correspond to 90% confidence level. 
"Fixed value indicates the Galactic absorption column density in units of 10^" cm^^. 
''Flux in units of 10""'^^ ergcm"^ s^^. 

230 PKS 0208-512 

231 The time averaged, background subtracted three XIS spectra of PKS 0208—512, when fitted 

232 jointly, are well described by a single absorbed power-law model, and the absorption column 

233 is consistent with the Galactic value A^h = 3.08 x lO^'^cm"^ (Dickey & Lockman 1990). 

234 We obtained the best fit photon index V = 1.68 ± 0.02 and the 2 — lOkeV fiux F2_iokeV = 

235 (1.37 ± 0.03) X 10~^^ergcm~^ s~^ with a chi-squared value of 0.91 for 250 dof. Figure[2] 

236 shows the spectra obtained with the XISs with residuals plotted against the best-fit power- 

237 law model with Galactic absorption. Although statistically acceptable, we notice that the 
23B residuals of the fits show moderate excess feature at low energies, below 1 keV. 

239 In the previous observation with BeppoSAX during a high fiux state (Tavecchio et 

240 al. 2002), -F2-iokcV ~ 4.7 x 10~^^ergcm^^ s^^ - which is a factor of three larger than in our 

241 Suzaku observations - the X-ray spectrum is well described by a power-law with photon index 

242 r ~ 1.7, similar to the Suzaku result. However, Tavecchio et al. reported that the spectrum 

243 was heavily absorbed below IkeV, indicating a column density of A^h = 1-67 x lO^^cm"^. 

244 FigureH] shows the Suzaku spectrum with residuals assuming such an increased value of A^h- 

245 The residuals indicate significant soft excess emission below 1 keV, if A^h is the same as found 

246 in the previous BeppoSAX observation. 

247 The variable soft X-ray emission of PKS 0208— 512 may indicate that the convex spec- 
24B trum observed by BeppoSAX refiects an intrinsic IC continuum shape, while the soft excess 

249 observed by Suzaku refiects the presence of an additional spectral component which becomes 

250 prominent when the source gets fainter (see Tavecchio et al. 2007; Kataoka et al. 2008). 
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Table 4: Results of the spectral fits to the Suzaku data with best fit models 



Object 


ModeP 




Phi^ 


-L low 


-F2-10koV'^ 


kT 
(keV) 


X? 


0208-512 


PL 


3.08 (fixed) 


1.68±0.03 




1.37±0.06 




0.91 (250) 


0827+243 


PL 


3.62 (fixed) 


1.46±0.04 




1.37±0.07 




0.84 (194) 


1127-145 


PL 


10.811J 


1.52±0.03 




3.36±0.08 




0.82 (330) 


1510-089 


PL+BB 


7.88 (fixed) 


1.32±0.03 




6.42±0.13 


0.15±0.03 


0.97 (405) 




PL+PL 


7.88 (fixed) 




2 85+0-88 






0.96 (405) 


3C 454.3 


PL 


-0.57 


1.62±0.01 




16.6±0.2 




0.97 (1089) 



"Spectral fitting models. PL, power-law function; PL+PL, double power-law function; PL+BB, power-law + 
blackbody model. 

''Differential spectral photon index. 

'^Differential spectral photon index at the low-energy X-ray band, when fitted with a double power-law func- 
tion. 

''Flux in units of 10^^^ ergcm^^ s^^. 

251 Therefore, to model in more detail the observed X-ray spectrum, we first considered a double 

252 power-law fit (PL -|- PL) in which the soft X-ray excess is represented by a steep power-law 

253 component. The absorption column is fixed at A^h = 1-67 x 10^^ cm~^, as given by Tavecchio 

254 et al. (2002). We obtained the photon indices Fi = 4.98 ± 0.30 and Fa = 1.7ll[]:[]i This 

255 provides an acceptable fit, with x^/dof = 0.90/248. We also considered an alternative fit 

256 consisting of a power-law function and a blackbody component. This model also gives a 

257 similarly good representation of the data with x^/dof = 0.90/248, implying F = 1.78 ± 0.03 

258 and the temperature of the introduced thermal component of kT = 0.092 ± 0.003 keV. Both 

259 fits appear to be as good as a single power-law with free absorption, and do not improve the 

260 goodness of fit. 

261 Q 0827+243 

262 The time averaged spectra of Q 0827-1-243 collected with the XISs are well fitted by an 

263 absorbed power- law model with a photon index F = 1.46 ±0.02 (x^ = 0.84 for 194 dof). The 

264 absorption column is consistent with the Galactic value of A^h = 3.62 x 10^° cm~^ (Dickey & 

265 Lockman 1990), and the flux over 2 - 10 keV is F2__iokcV = (1.37±0.04) x 10"^^ g-i^ 

266 As shown in Figure[2|, there is no evidence for any additional spectral feature in the soft band. 

267 This result is in good agreement with previous Chandra observations of the core (Jorstad & 
26B Marscher 2004), revealing that the X-ray continuum is well described by a power-law model 
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269 (r ~ 1.4) with Galactic absorption. 

270 PKS 1127-145 

271 We first fitted the XISs spectra with a single power-law model with a Galactic absorption of 

272 Nil = 3.83X 10^0 cm-2 (Murphy et al. 1996). We obtained the photon index of T = 1.41±0.01 

273 Xr = 1-01 331 dof), but the residuals show a substantial deficit of photons at low energies 

274 (FigureEI). To investigate this deficit in more detail, we fitted the spectra with a single 

275 power-law and a free absorption model. This model represents well the spectra with the best 

276 chi-squared value of 0.82 for 330 dof (Figure|3]), indicating that the column density is higher 

277 than the Galactic value at the 99.9% confidence level. For this model the photon index is 
27B r = 1.51±0.02 and the unabsorbed X-ray flux is -F2„iokcV = (3.36±0.05) x 10~^^ ergcm"^ s~^ 

279 The best-fit column density is = (1.08 ± 0.09) x 10^^ cm~^, which is similar to the one 

280 found in previous Chandra and XM.M- Newton observations (A^h ~ 1-2 x 10^^ cm~^) during a 

281 high state with F2_iokeV ~ 6 x 10~^^ ergcm"^ s"^ (Bechtold et al. 2001; Foschini et al. 2006). 

282 We note that the Galactic absorption and a broken power-law model also well represents the 

283 spectra with x^/dof of 0.81/329. In this model, the spectrum below £^brk = 1-50 ± 0.08 keV 

284 is rather hard (Fi = 1.10 ± 0.05), and the high energy photon index is r2 = 1.50 ± 0.02. 

285 PKS151 0-089 

286 Figure^] shows the XISs and HXD/PIN spectra of PKS 1510-089 (including residuals), plot- 

287 ted against the best-fit power-law model with Galactic absorption, using the overall X-ray 

288 data between 0.3 and 40keV. The best fit photon index is F = 1.37 ± 0.01 and the unab- 

289 sorbed X-ray flux is F2_iokeV = (6.31 ± 0.07) x 10~^^ergcm~^ s^^ However, this model did 

290 not represent the spectra well yielding a chi-squared value of 1.06 for 407 dof. The residuals 

291 indicate some excess emission at low energies. 

292 To represent the observed X-ray spectra, we tried the same analysis as for PKS 0208—512. 

293 We first fitted the data by a double power-law model with Galactic absorption. We obtained 

294 the photon indices Fi = 2.SA'^^^ and F2 = 1.26;!;o;Qg. This provides an acceptable fit, with 

295 x?/dof = 0.96/405. The improvement of the chi-squared statistic is significant at more than 

296 the 99.9% confidence level when compared to the single power-law model. Next, we consid- 

297 ered an alternative fit consisting of a power-law function and a blackbody component. This 

298 model also gives a good representation of the data, with of 0.97 for 405 dof, indicating that 

299 the photon index is F = 1.32 ± 0.02 and the temperature of the introduced thermal compo- 

300 nent is kT = 0.15±0.02keV. This result is consistent with previous Suzaku (F = 1.24±0.01; 

301 Kataoka et al. 2008) and BeppoSAX observations (F = 1.39 ± 0.08; Tavecchio et al. 2000). 
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302 3C 454-3 

303 We first fitted the XISs and PIN spectra with a single power-law model with a Galactic 

304 absorption of A^h = 7.24 x 10^° cm~^ (Murphy et al. 1996). We obtained the photon index of 

305 r = 1.41 ±0.01 {xf = 1-00 1090 dof), but the residuals show some scatter around 1 keV. 

306 To investigate this scatter in more detail, we fitted the spectra with a single power-law and 

307 a free absorption model. This model represents well the spectra with the best chi-squared 

308 value of 0.97 for 1089 dof, indicating that the column density is higher than the Galactic 

309 value at the 99.9% confidence level. For this model the photon index is F = 1.51 ± 0.02 and 

310 the unabsorbed X-ray flux is -F2-iokcV = (3.36 ± 0.05) x lO^^^ergcm"^ s~^. 

311 In the case of 3C 454.3 we obtained the best flt to the XISs and HXD/PIN spectra 

312 assuming an absorbed power-law model with a photon index of F = 1.61 ± 0.01 and a 

313 column density of A^h = (9.07 ± 0.35) x lO^'^cm"^, which is larger than the Galactic value 

314 at the 99.9% confldence level. The unabsorbed 2 — lOkeV flux is -F2_iokov = (1-66 ± 0.01) x 

315 10~^^ ergcm"^ s~^ (Figure [2]). The spectra can be fltted with both the Galactic absorption 

316 and a broken power-law model as well as the above model {x^/dof of 0.96/1088). In the 

317 former case, the photon indices are Fi = 1.471q'q3 ^2 = 1.61 ± 0.01, while the break 
31B energy is i^brk = 1-291q°^ keV. In addition, we reanalyzed the previous Suzaku data collected 

319 in December 2007 during the high state (Donnarumma et al. 2010). The time averaged 

320 XISs and HXD/PIN spectra was well described by a single absorbed power-law model with 

321 F = 1.64 ± 0.01, implying the flux F2_iokov = (3.09 ± 0.02) x 10"^^ ergcm'^ s-\ which is 

322 larger by a factor of two than the one found in our 2008 observations. The absorption column 

323 also shows a higher value of A^h = (1.07 ± 0.03) x 10^^ cm~^. 

324 Figure|5] shows the unfolded spectra obtained in 2007 (high state) and 2008 (this work; 

325 low state). The bottom panel shows the residuals by subtracting the spectra in the high 

326 state from those in the low state. The excess emission at low energies is clearly visible in 

327 the residuals. 

32B The previous X-ray observations of 3C 454.3 often indicated some additional absorption 

329 in excess to the Galactic value. For example, Villata et al. (2006) reported A^h = (1.34 ± 

330 0.05) X 10^^ cm~^ in the Chandra data collected in May 2005, during the outburst phase 

331 (-F2-8keV ~ 8.4 X 10~^^ ergcm"^ s~^, which is ~ 5 times higher than in our observation). An 

332 even higher hydrogen column density was found by Giommi et al. (2006), when fltting the 

333 April-May 2005 data taken by the Swift XRT (A^h ~ 2 - 3 x 10^^ cm'^), and by Atari et al. 

334 (2007, 2008), using the July and December 2006, and May 2007 data taken by XyiM-Newton. 

335 Assuming that the intrinsic absorption in 3C 454.3 is the same as reported in Villata et al. 

336 (2006), we flt our Suzaku data flrst by a double power-law function, obtaining Fi = 3.661q 33 

337 and F2 = 1.61 ± 0.02. This provides an acceptable flt, with x^/dof = 0.97/1088. Next 
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33B we consider an alternative fit consisting of a power-law and a blackbody component. This 

339 model gives a good representation of the data, with x^/dof = 1.00/1088, a photon index 

340 of r = 1.65 ± 0.01, and a temperature of kT = 0.105 it 0.004keV. However, the fitting 

341 results do not improve the goodness of fit compared with the single power-law model with 

342 free absorption. 

343 3.1.3. Time-resolved spectral analysis 

344 In order to investigate the X-ray spectral evolution of each object, we divided the total 

345 exposure into one-orbit intervals (~ 5760 s). We fitted the overall XIS spectra between 0.3 

346 and 10 keV with an absorbed simple power-law function. The photoelectric absorbing column 

347 densities were fixed at the values derived in §3.1.2. Figure[6] shows the relation between 

348 the 2 — 10 keV fluxes versus the photon indices measured by the Suzaku XISs. Significant 

349 spectral variation is seen in PKS 0208-512 (F = 1.4-1.8), Q 0827+243 (F = 1.2-1.6), 

350 PKS 1127-145 (r = 1.4 - 1.6), and PKS 1510-089 (F = 1.3 - 1.5). In the case of 3C 454.3, 

351 the X-ray photon index is only weakly variable around the mean value F ~ 1.6. Figure[6] 

352 clearly reveals a spectral evolution with the X-ray spectra hardening as the sources become 

353 brighter. Such a trend is often observed in high-frequency-peaked BL Lac objects (e.g., 

354 Kataoka et al. 1999), but it has never been observed so clearly in FSRQs (but see Kataoka 

355 et al. 2008 for PKS 1510-089). 

356 3.2. Swift 

357 Since the effective area of the Swift XRT is less than 10% of the Suzaku XIS in the 

358 0.5—10 keV range, detailed spectral modeling is difficult using Swift data. Furthermore, the 

359 average exposure for the Swift observation was only a few kiloseconds, which was much less 

360 than the Suzaku exposure. We therefore fit the XRT data simply with a power-law model 

361 with Galactic absorption in the energy range 0.3—10 keV for the cross-calibration between 

362 the two instruments. The results of the spectral fits are summarized in Tabled We can see 

363 that the results obtained with Suzaku and Swift are consistent within the range of error 

364 except PKS 1127-145. 

365 The UVOT fluxes in each filter were corrected for Galactic extinction following the 

366 procedure described in Cardelli et al. (1989). We generated a list of the amount of extinc- 

367 tion that needs to be accounted for in each filter. Ax = Es-vio-Rv + b), where a and b 

368 are constants. The Cardelli procedure provides a good approximation to the UV-through- 
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Table 5: Results of the spectral fits to the Swift XRT data using a power-law with Galactic 
absorption 



Object 



-r2-10kcV Xi 



(dof) 



0208-512 
0827+243 
1127-145 
1510-089 
3C 454.3 



3.08 (fixed) 
3.62 (fixed) 
3.83 (fixed) 
7.88 (fixed) 
7.24 (fixed) 



1.96±0.24 
1.46±0.35 
1.28±0.03 
1.38±0.08 
1.53±0.09 



1.37±0.06 

1 1 O+0.55 

6.14±0.23 

"•^^-0.62 

17.6±2.1 



0.33 (13) 
1.12 (8) 
1.28 (94) 
0.71 (31) 
1.41 (18) 



Errors correspond to la confidence level. 
"Fixed value indicates the Galactic absorption column density in units of 10^" cm^^. 
''Flux in units of 10""'^^ ergcm"^ s^^. 

369 IR Galactic dust extinction as a function of the total-to-selective extinction, Ry, which 

370 throughout this paper we assume to be Ry = 3.1, which is the mean Galactic value. The 

371 observed magnitudes and correction factors for each of the filters are summarized in Table[6] 
and TableO respectively. 

Table 6: Swift UVOT magnitudes of five FSRQs 



Object 


V 


b 


u 


uvwl 


uvm2 


uvw2 


E(B-V) 




(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 




0208-512 


17 64+"- 1^ 




1 7 1 7+0.09 
^'■^'-0.08 


16.84+0.07 


16.71+0.07 


17.001-.- 


0.022 


0827+243 






16.57+0.01 








0.033 


1127-145 


16.48+0.02 


16.70+0.01 


15.64+0.01 


15.51+0.01 


15.56+0.01 


15.79+0.01 


0.037 


1510-089 






16.01+0.02 


16.27+0.02 


16.13+0.02 




0.097 


3C 454.3 


16 05+°-°^ 

iO.UO_g og 


16.52+0.06 


15.72+0.06 


15.75+0.06 


15.81+0.08 


16.06+0.05 


0.107 



Observed magnitude for each observation using specific filter (Galactic extinction not corrected). 



3.3. Fermi LAT 



374 To study the average spectra of five objects during the four or five months of observa- 

375 tions, we use the standard maximum-likelihood spectral estimator provided with the LAT 
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Table 7: Correction factors for the Galactic extinction in UV and optical filters 



Par am 


V 


b 


u 


uvwl 


uvni2 


uvw2 




(nm) 


547 


439 


346 


260 


249 


193 






1.0015 


0.9994 


0.9226 


0.4346 


0.3494 


-0.0581 






0.0126 


1.0171 


2.1019 


5.3286 


6.1427 


8.4402 




0208-512 


0.07 


0.09 


0.11 


0.15 


0.16 


0.18 




0827+243 






0.16 








Ax"" 


1127-145 


0.12 


0.15 


0.18 


0.25 


0.27 


0.31 




1510-089 






0.48 


0.65 


0.70 






3C 454.3 


0.33 


0.44 


0.53 


0.71 


0.77 


0.88 



"Center wavelength for each optical and UV filter. 

^Parameters for calculating Galactic extinction for optical and UV filters, calculated according to the pre- 
scription in Cardelli et al. (1989). The Galactic reddening was taken from Schlegel et al. (1998). 

376 science tools gtlike. This fits the data to a source model, along with models for the uniform 

377 extragalactic and structured Galactic backgrounds. Photons were extracted from a region 

378 with a 10° radius centered on the coordinates of the position of each object. The Galactic 

379 diffuse background model is the currently recommended version (glLiem_v02 0), with the 

380 normalization free to vary in the fit. The response function used is P6_V3_DIFFUSE. 

381 For simplicity, we model the continuum emission from each source with a single power- 

382 law. It is likely that such a model might be too simple, as shown in the paper reporting 

383 spectra of bright Fermi blazars (Abdo et al. 2010), where the gamma-ray data suggest a 

384 steepening of the spectrum with energy, well- described as a broken power-law. However, here, 

385 we are reporting cases of blazars in low- level activity states and thus relatively faint, where 

386 fits to a broken power-law model would result in poorly constrained spectral parameters for a 

387 more complex model; furthermore, we note that the use of such more complex spectral model 

388 in the gamma-ray band does not alter our conclusions or significantly change the parameters 

389 in Table 10. The extragalactic background is assumed to have a power-law spectrum, with 

390 its spectral index and the normalization free to vary in the fit. From an unbinned gtlike 

391 fit the best fit photon indices are F = 2.33 ± 0.05 for PKS 0208-512, F = 2.62 ± 0.35 for 

392 Q 0827+243, F = 2.77 ± 0.14 for PKS 1127-145, F = 2.48 ± 0.03 for PKS 1510-089, and 

393 F = 2.51 ± 0.02 for 3C 454.3 (see also Table[8]). Here only statistical errors are taken into 



^This model is available for download from the Fermi Science Support Center, 
http://fermi.gsfc.nasa.gov/ssc 
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394 account, and we report fluxes using spectra extrapolated down to 100 MeV. In the case of 

395 bright sources (PKS 1510—089 and 3C 454.3), we also analyzed the data collected during the 

396 Suzaku observing period to construct the simultaneous broad-band spectra spectra. 

397 Table [9] summarizes the flux in seven energy bands obtained by separately running 
39B GTLIKE for each energy band; 200-400 keV, 400-800 keV, 800-1600 keV, 1600-3200 keV, 
399 3200-6400 keV, 6400-12800 keV, 12800-25600 keV, respectively. 

Table 8: Results of the spectral flts to the Fermi LAT data 



Object r i^>100Mev" TS 



0208-512 2.33±0.05 0.26±0.03 1484 

0827+243 2.62±0.36 0.05±0.04 58 

1127-145 2.75±0.14 0.15±0.04 234 

1510-089 2.48±0.03 0.69±0.04 4224 

3C 454.3 2.50±0.02 2.55±0.08 25144 

1510-089^= 2.28±0.27 0.91±0.51 59~" 

30454.3'= 2.62±0.13 2.59±0.58 281 



"Flux in units of 10 ^phcm ^ s ^. 

^Test statistic: defined as TS = 2(logL — logLo)i where L and Lq are the Ukehhood when the source is 
included or not. 

^Oorresponding data collected during the Suzaku observing period. 



400 4. Discussion 

401 4.1. Broad-Band Spectra spectral fits 

402 We constructed the broad-band spectral energy distribution (SED) ranging from the 

403 radio to 7-ray bands for the flve observed FSRQs, and these are shown in FigureH Here 

404 the fllled red circles and solid lines represent simultaneous data from the UV/optical {Swift 

405 UVOT), X-ray (Suzaku) and 7-ray [Fermi LAT) observations. Quasi- simultaneous data 

406 are also shown as red open triangles and dashed lines. Historical radio (NED) and 7- 

407 ray (EGRET) data are also plotted as fllled blue circles. Green symbols in the SEDs of 

408 PKS 1510—089 and 3C 454.3 denote the previous simultaneous observations (Kataoka et al. 

409 2008; Donnarumma et al. 2010). 

410 In order to model the constructed SEDs, we applied the synchrotron-inverse Compton 
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Table 9: Results of Fermi LAT data analysis from 200 MeV to 25600 MeV (Flux in units of 
lO-^phMeV-^cm-^s-i) 



Object 


Band I'' 


Band 2^ 


Band 3^ 


Band 4^ 


Band 5"" 


Band 6^ 


Band 7^ 


0208-512 
0827+243 
1127-145 
1510-089 
3C 454.3 


3.24±0.25 
0.73±0.22 
2.65±0.34 
9.74±0.38 
37.22±0.67 


3.69±0.30 
1.03±0.28 
2.50±0.43 
10.62±0.48 
39.36±0.91 


4.19±0.42 

1.08±0.36 

3.21±0.68 

10.19±0.67 

43.79±1.42 


3.81±0.61 

0.92±0.50 

3.17±1.05 

11.57±1.09 

44.20±2.31 


2.88±0.79 

7.51±1.47 
26.65±2.96 


1.95±1.00 

7.15±2.26 
17.98±3.99 


1.51±1.49 

5.14±3.01 
4.55±3.28 



"Band 1: 200-400 MeV, Band 2: 400-800 MeV, Band 3: 800-1600 MeV, Band 4: 1600-3200 MeV, Band 
5: 3200-6400 MeV, Band 6: 6400-12800 MeV, and Band 7: 12800-25600 MeV 

411 (IC) emission model described in Tavecchio & Ghisellini (2008), where both synchrotron and 

412 external (BLR and DT) photons are considered as seed radiation fields contributing to the IC 

413 process (SSC+ECR). The electron distribution is modeled as a smoothly broken power-law: 

(\ ni-n2 
1 + — , (1) 

7br/ 

414 where K (cm~^) is a normalization factor, ni and n2 are the energy indices below and 

415 above the break Lorentz factor 7br. The electron distribution extends within the limits 
«6 7min < 7 < 7max- Wc also assume that the 'blazar emission zone', with the comoving size R 

417 and magnetic field intensity B, is located at the distance r such that vq < r < tblr < ''"dt, 

418 where tq is the distance below which the photon energy density in the jet rest frame is 

419 dominated by the direct radiation of the accretion disk, tblr is the characteristic scale of 

420 the broad-line-region, and tdt is the scale of the dusty torus (see the discussion in Tavecchio 

421 & Ghisellini 2008 as well as in Sikora et al. 2009). This choice, while somewhat arbitrary, 

422 has been validated by a number of authors modeling broad-band spectra of FSRQs. Hence, 

423 the comoving energy density of the dominant photon field — provided by the BLR — is 



424 where Tj is the jet bulk Lorentz factor, and the BLR is assumed to reprocess t^blr — 10% 

425 of the disk luminosity Lj. Finally, we assume that the jet viewing angle is in all the cases 

426 6j ^ l/^j, so that the jet Doppler factor 6j ~ Tj. 

427 The results of model fitting are shown in different panels of Figure[71 and the resulting 
42B parameters are summarized in Table[TOl In the context of this model, where we assume the 
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429 dissipation region to be between the immediate vicinity of the accretion disk but within the 

430 BLR, it is clear that in all cases the LAT fluxes are dominated by the IC/BLR component, 

431 while in the X-ray band both IC/BLR and IC/DT processes may contribute at a comparable 

432 level. In addition, the SSC emission seems negligible, being in particular too weak to ac- 

433 count for the soft X-ray excess discussed in the previous sections. This excess, on the other 

434 hand, may be well represented by the high-energy tail of the synchrotron continuum, or an 

435 additional blackbody-type spectral component. 

Table 10: Model parameters used to calculate the SEDs of five FSRQs 



Object 


ni 


n2 


Tmiii 


7br 


Tmax 


K 


r, 


R 


B 




fBLR 


r-DT 














[lO^cm-^] 




[10^^ cm] 


[G] 


[10^6 erg/s] 


[IQi^cm] 


[10^8 cm] 


0208-512 


2 


3.3 


3.0 


700 


4.3 xlO'^ 


2.2 


15 


1.8 


1.1 


1.5 


0.76 


3.0 


0827+243 


2 


3.3 


1.5 


300 


1.0 xlO" 


8.5 


10 


1.8 


3.8 


2.0 


1.3 


4.2 


1127-145 


2 


3.4 


1.2 


110 


5.0 xlO^ 


0.65 


10 


6.5 


4.1 


10 


0.8 


10 


1510-089 


2 


3.5 


3.0 


190 


4.6 X 10-* 


0.81 


13 


4 


0.8 


0.3 


0.48 


4 


3C 454.3 


2 


3.8 


1.0 


290 


3.0 X 10'' 


4.5 


12 


3.2 


0.8 


4.0 


1.5 


30 



436 Based on the model results, for each object we compute the ratio of the comoving energy 

437 densities stored in jet electrons and the magnetic field, 

43B where B is the magnetic field intensity in the emission region. In addition, we compute the 

439 implied total kinetic jet power as 

L, = T^R'cT] + U's + U;) , (4) 

440 where R is the emission region linear size, and is the energy density of cold protons. 

441 The latter parameter is estimated assuming one proton per ten electron-positron pairs (see 

442 the discussion in Sikora et al. 2009), namely U' = O.lmpC^ P'^'^'^ N'H) d'j. The resulting 

443 total kinetic power of the outflow is then compared with the accretion luminosity (assuming 

444 standard accretion disk with 10% radiative efficiency), by means of the evaluated efficiency 

445 parameter rjj = Lj/L^cc — Lj/lOLf^, where is the disk luminosity implied by the model 

446 fitting (see TablefTOj) . Note that with the above model assumptions and the model parameters 

447 inferred by us, the jets of objects considered here are dynamically dominated by cold protons, 
44B Up/U'g ~ 200/(7) > since the mean Lorentz factor of the radiating ultra-relativistic 
449 electrons is in all the cases (7) ^ 200 (see TablefTT]) . 



(3) 
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Table 11: Jet parameters of five FSRQs in low-activity states 



Object 




[lO^^ erg/s] 




(7) 


0208-512 


2 


0.8 


0.06 


16 


0827+243 


0.6 


2.8 


0.14 


8 


1127-145 


0.03 


5.8 


0.06 


5 


1510-089 


0.9 


1 


0.35 


12 


3C 454.3 


7 


9.4 


0.23 


5 



450 Some of the derived jet parameters for five luminous blazars in their low-activity states 

451 are significantly different from the analogous parameters claimed for the flaring states, even 

452 in the same object. For example, in the case of the high-activity state of PKS 1510—089, 

453 Kataoka et al. (2008) estimated (under the same assumptions regarding the jet content as in 

454 this paper) the total kinetic power of the jet as Lj ~ 2.7 x 10'^^ ergs^^, which is larger than 

455 the value derived in this paper, by about a factor of 3. In addition, our model values of the 

456 jet bulk Lorentz factors are also systematically lower than the ones given in the literature 

457 (Tj ~ 10 versus 20). Interestingly, other jet parameters, such as magnetic field intensity 

458 - S ~ IG - and the equipartition ratio, U'^/U'^ ~ 1, or the general spectral shape of the 

459 electron energy distribution, are comparable to the ones found for fiaring FSRQs, (albeit 

460 with a substantial scatter). It should be noted in this context, however, that for the three 

461 sources considered in this paper (namely PKS 1127—145, PKS 1510—089, and 3C 454.3), the 

462 flaring states were analyzed in a framework of the IC/DT model (Blazejowski et al. 2004, 

463 Kataoka et al. 2008, Sikora et al. 2008, respectively), while here, we argue that the IC/BLR 

464 contribution is dominant, as motivated by the detected relatively short (day) variability 

465 timescale of the X-ray continua. On the other hand, as discussed recently in Sikora et al. 

466 (2009), there is so-called a 'conspiracy' between the IC/BLR and IC/DT models, in a sense 

467 that the resulting inferred jet parameters are comparable in both cases. Hence, we can 

468 safely conclude that the low- and high-activity states of luminous blazar sources are due 

469 to the low and high total kinetic power of the jet, respectively, possibly related to varying 

470 bulk Lorentz factors within the blazar emission zone. And indeed, keeping in mind that the 

471 highly dynamical and complex jet formation processes in the closest vicinity of supermassive 

472 black holes - most likely shaped by accretion process subjected to several possible instability 

473 of the jet fuel, especially when the accretion rate is close to Eddington - such a signiflcant 

474 variation in the total kinetic output of the outflow should not be surprising. Further support 

475 for this scenario comes from the fact that the jet efficiency factors estimated here, r]j < 1, 

476 are significantly lower than the ones found for powerful blazars in their flaring states (see 

477 Sambruna et al. 2006, Ghisellini et al. 2009), even if the difference in the jet proton content 
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478 adopted by various authors is taken into account. 



4.2. Spectral Evolution 



480 
481 



As shown in § 3.1.3, the X-ray spectra of the FSRQs analyzed here flatten with increasing 
flux. For Fj ~ 5j ~ 10 and the dominant IC/BLR emission process, the electrons emitting 

482 the observed 1 — 10 keV photons have Lorentz factor 7 ~ 7min ~ few. The electrons emitting 

483 X-ray photons in these sources are very low-energy, so cooling effects cannot play any role 

484 in the observed spectral evolution. In particular, it can be easily demonstrated that in a 

485 framework of our model (i.e., for the dominant IC/BLR energy losses), a strong cooling 

486 regime is expected only for the electrons with Lorentz factors greater than 



496 
497 
498 



501 



7c 



37r meC^ r^LR 



O-T R F j T^BLR^d 

^ 350 r ''^^^ V P^^^rT^ V7 VV ^ V (5) 

VlOiScm/ V 0.1 / \10 J VlO^^erg/sy VlO^^cmy ' 

487 This, for the fitting parameters as given in TablelTOl is typically above or just around the 

488 break Lorentz factor, 'jcr ^ 7br (in agreement with the discussion in Sikora et al. 2009). 

489 Adiabatic losses, if present, should not result in changing the slope of the power-law X-ray 

490 continua as well. Thus, one may suspect that the revealed spectral changes are shaped 

491 by the acceleration process within the blazar emission zone. In the case of relativistic jets 

492 the relevant acceleration processes are still quite uncertain, although, as pointed out by 

493 Kataoka et al. (2008) and Sikora et al (2009), the repeatedly observed flat X-ray photon 

494 indices F < 1.5 seem to favor the mechanism discussed by Hoshino et al. (1992) for the low- 

495 energy segment of the electron energy distribution. In this model, the low-energy electrons 
(with Lorentz factors, roughly, 7 < mp/rrie) are accelerated by a resonant absorption of the 
cyclotron emission generated by cold protons reflected from the shock front. As shown later 
by Amato & Arons (2006), the power-law slope of these accelerated electrons depends on 

499 the relative number of electrons to protons at the shock front. Hence, a larger fraction of 

500 the energy carried by jet protons during the higher-activity states should in principle result 
in a more efficient acceleration of jet electrons and their flatter spectrum, in agreement with 



502 the observed X-ray spectral evolution discussed here. 

503 The above interpretation, on the other hand, would imply a signiflcant variability in the 

504 7-ray frequency range. Indeed, the broken power-law form of the electron energy distribution 

505 revealed by our spectral modeling discussed in the previous section implies the 7-ray flux 
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506 = [vFy]^ around the IC spectral peak z/^ ~ 10^^ Hz should be, roughly 

F.^Fx ':^10^(^-^)Fx, (6) 

507 where Fx is the monochromatic X-ray flux measured around vx ~ 10^^ Hz, and F is the 

508 observed X-ray photon index. For example, our analysis for PKS 0208— 512 indicates a 

509 photon index Fi ~ 1.8 for an X-ray flux Fx,i ~ 1-2 x 10~^^ergcm~^ s~^ in the lower state, 

510 and F2 ~ 1.5 for Fx,2 ~ 1-6 x 10~^^ ergcm"^ s~^ in the higher state. Thus, if the observed X- 

511 ray variability is due to flattening of the electron energy distribution during the acceleration 

512 process, one should observe the 7-ray variability of the order of 

:^^10^(ri-r,):^_20. (7) 
F^, 1 Fx, 1 

513 However, during the simultaneous Fermi observation, no significant 7-ray variability was 

514 observed for the analyzed sources, at least within one day timescale. 

515 Therefore, the most viable explanation for the observed X-ray spectral evolution is that 

516 the IC power-law slope remains roughly constant during the flux variations, but the amount 

517 of contamination from the additional soft X-ray component increases at low flux levels, 

518 affecting the spectral fitting parameters at higher photon energies (> 2keV). Note that in 

519 such a case the expected gamma-ray variability should be of the same order as the X-ray 

520 variability, namely F^, 2/-P7, 1 — Fx, 2/ Fx, i ~ 1-3. 

521 We finally note in this context that, as shown in §3.1.2, the previous BeppoSAX data 

522 for PKS 0208— 512 collected during the high state indicated a convex X-ray spectrum, and 

523 an excess absorption below 1 keV with a column density of A^h ~ 1-67 x 10^^ cm"^ exceeding 

524 the Galactic value by more than a factor of 5. However, the X-ray photon index was similar 

525 to the one implied by our Suzaku observations (F ~ 1.7). Therefore, the convex spectrum 

526 observed by BeppoSAX may reflect an intrinsic shape of the IC emission involving the low- 

527 energy cut-off in the electron energy distribution around 7 ~ 1, as expected in the EC/BLR 

528 model (Tavecchio et al. 2007), which is only diluted during the low-activity states due to 

529 the presence of an additional soft X-ray spectral component. 

530 Similar trend has been observed in 3C 454.3. To illustrate this, in FigureElwe selected 

531 the data which have a similar power-law slope (F ~ 1.6) and plotted the absorption column 

532 versus 2 — 10 keV flux densities derived from the Chandra (Villata et al. 2006), Swift (Giommi 

533 et al. 2006), XMM- Newton (Raiteri et al. 2007, 2008), and Suzaku (this work) observations. 

534 We can see that there is a trend of increasing the absorption value with source brightness, as 

535 previously reported by Raiteri et al. (2007; 2008). These results may again be explained by 
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536 the soft excess emission being more important when the source gets fainter, and becoming 

537 almost completely "hidden" behind the hard X-ray power-law when the source gets brighter. 

53B From the spectral fitting of the Suzaku data, we showed in §3.1.2 that the soft X-ray 

539 excess may be represented either by a steep power-law (F ~ 3 — 5) or a blackbody-type 

540 emission {kT ~ 0.1 — 0.2 keV). Since the synchrotron peak of each source is located around 

541 optical photon energies (see Figure[7]), the high-energy synchrotron tail may possibly account 

542 for the observed soft X-ray excess emission, especially if being modified by the Klein- Nishina 

543 effects (see the discussion in Sikora et al. 2009 and Kataoka et al. 2008). On the other hand, 

544 the bulk-Compton spectral component produced by Comptonization of the UV accretion 

545 disk by cold electrons in the innermost parts of relativistic jets (e.g., Begelman & Sikora 

546 1987) is a natural explanation for the apparent soft X-ray excess component. 

547 5. Conclusions 

548 We have presented the observations and analysis of the data for the 7-ray-loud blazars, 

549 PKS 0208-512, Q 0827+243, PKS 1127-145, PKS 1510-089, and 3C 454.3, obtained with 

550 the Suzaku, Swift UVOT and Fermi LAT. Observations were conducted between 2008 Octo- 

551 ber and 2009 January. These observations allowed us to construct broadband spectra of the 

552 sources in the low 7-ray activity state, covering optical to GeV photon energy range. Our 

553 results follows: 



554 1. The X-ray spectra of five FSRQs are well represented by an absorbed hard power-law 

555 model (F ~ 1.4 - 1.7). For PKS0208-512, PKS 1127-145, and 3C454.3, the fitted 

556 absorption column is larger than the Galactic value (but we note that the "excess 

557 absorption" is not a unique representation of X-ray spectra of those blazars). Com- 

558 pared with previous X-ray observations, we see a trend of increasing apparent X-ray 

559 absorption column with increasing high-energy luminosity of the source. 

560 2. Suzaku observations reveal spectral evolution of the X-ray emission: the X-ray spec- 

561 trum becomes harder as the source gets brighter. Such spectral changes are most 

562 likely due to the underlying and steady low-energy spectral component which becomes 

563 prominent when the inverse-Compton emission gets fainter. This soft X-ray excess can 

564 be explained as a contribution of the high-energy tail of the synchrotron component, 

565 or bulk-Compton radiation. 

566 3. We adopt the location of the blazar emission region to be outside of the immediate 

567 vicinity of the accretion disk but within the BLR, and within the context of this 
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568 model, we find that the contribution of the synchrotron self-Compton process to the 

569 high-energy radiative output of FSRQs is neghgible even in their low-activity states. 

570 4. We argue that the difference between the low- and high-activity states in luminous 

571 blazars is due to the different total kinetic power of the jet, most hkely related to 

572 varying bulk Lorentz factor of the outflow within the blazar emission zone. 
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Fig. 1. — Light curves of five FSRQs: 0.5 — 2 keV [upper panels), 2 — 10 keV [middle panels), 
and 0.5 — lOkeV [bottom panels). All the light curves were binned at 5760s, corresponding 
to the period of the Suzaku orbit. 
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Fig. 2. — Suzaku spectra of five FSRQs: the top panel shows the data, plotted against a 
power-law model with the Galactic absorption. The bottom panel shows the residuals for 
the power-law fit. For 0208—512 and 1510—089, the data below 1 keV are in excess to the 
model. On the other hand, for 1127—145, the residuals show a substantial deficit of photons 
at low energies. For 3C454.3, some scatter around 1 keV in the residual panel is seen. 
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Fig. 3.— Best fit Suzaku spectrum of PKS 1127-145, PKS 1510-089 and 3C454.3. The 
top panel shows the data, plotted against an absorbed power-law model. The bottom panel 
shows the residuals to the power-law fit. 
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Fig. 4. — Suzaku spectrum of PKS 0208—512, with residuals assuming a column density of 
Nil ~ 1-67 X 10^^ cm^^. Deviations due to soft excess emission can be clearly seen. 
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Fig. 5. — Unfolded spectra of 3C 454.3 obtained in 2007 (high) and 2008 (low), respectively. 
The bottom panel shows the residuals by subtracting the spectra in the high state from those 
in the low state. 
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Fig. 6. — Correlation of the 2 — lOkeV flux vs. photon index of five blazars as measured by 
the Suzaku XISs. 
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Fig. 7. — Overall SED of five sources constructed with broad-band data obtained during 
2008 October to 2009 January (filled red circles and solid bow-tie). Quasi simultaneous data 
are also shown (open red circles). Historical radio (NED) and 7-ray (EGRET) data are also 
plotted as filled blue circles. Green symbols in the SEDs of PKS 1510— 089 and 3C 454.3 
denote the previous simultaneous observations (Kataoka et al. 2008; Donnarumma et al. 
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Fig. 8. — Fluxes in the 2 — lOkeV band for different observations of 3C 454.3 vs. A^h for the 
fits with an absorbed power-law model. The dashed line indicates the Galactic absorption 
column. This figure indicates that the intrinsic X-ray spectrum is not a simple power-law, 
but instead, it shows some curvature, which may depend on the X-ray brightness. 



